Magnetization Jump in a Model for Flux Lattice Melting at Low 

Magnetic Fields 



Seungoh Ryu 
D. Stroud 

Department of Physics, Ohio State University, Columbus, OH 43210 

(February 1, 2008) 

Abstract 

Using a frustrated XY model on a lattice with open boundary conditions, 
we numerically study the magnetization change near a flux lattice melting 
transition at low fields. In both two and three dimensions, we find that the 
melting transition is followed at a higher temperature by the onset of large 
dissipation associated with the zero-field XY transition. It is characterized by 
the proliferation of vortex-antivortex pairs (in 2D) or vortex loops (in 3D). At 
the upper transition, there is a sharp increase in magnetization, in qualitative 
agreement with recent local Hall probe experiments. 
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Among the most debated aspects of high-T c superconductors is the nature of the su- 
perconducting transition in a magnetic field. It is widely believed that in the presence of 
disorder, the transition may fall into the universality class of a vortex glass JI|, or a Bose 
glass 0. In addition, disorder may produce an entirely new transition characterized by the 
disappearance of the Bragg peaks of the Abrikosov lattice above a certain magnetic field 
||. But even in the absence of disorder, the phase transition is not fully understood. At 
high applied fields H parallel to c, several numerical studies strongly suggest that there is 
a first-order phase transition with a finite latent heat |§-§|], and a corresponding magneti- 
zation discontinuity ||. At lower magnetic fields, the evidence is more ambiguous. At such 
fields, the thermodynamics of the field-induced vortices should be strongly influenced by 
fluctuations of the zero-field XY degrees of freedom, which are far more numerous than the 
field induced vortex degrees of freedom |7] . 

Experimental evidence for a first-order transition with a magnetization jump has been 
reported for both Bi 2 Sr 2 CaCu208+5 || and YBa 2 Cu30 7 _ ( 5 ||. These measurements show 
that the transition is accompanied by an increase in vortex number density, suggesting a 
melting analogous to that of water. 

In this paper, we present model calculations of the superconducting-normal (SN) transi- 
tion in both two and three dimensions (2D and 3D) in the absence of disorder. Our results 
strongly suggest that in both cases, the SN transition at low B occurs in two distinct stages. 
The first is the melting of the vortex lattice at T m , and is characterized by disappearance of 
the Bragg peaks associated with solid order. The second stage is driven by the XY degrees 
of freedom, and occurs at temperatures above T m . That zero-field XY transition is accom- 
panied by the unbinding of vortex-antivortex pairs (in 2D), or proliferation of vortex loops 
(in 3D) ]HJ. In either case, these excitations screen the repulsive interactions between the 
field-induced vortices or vortex lines |11|] . Consequently, the density of field-induced vortices 
becomes larger, leading to an increase in the absolute magnetization as observed experi- 
mentally. In 3D, this screening is accompanied by the vanishing of 7 2Z , the z-component of 
the helicity modulus, which measures the rigidity of the phase degrees of freedom in the c 

2 



direction |JL2]| . Since this occurs at a sharp temperature, the second stage could be a true 
phase transition in 3D. 

Our calculations are based on the widely-studied frustrated XY model on a stacked 
triangular lattice 0. This model has the Hamiltonian TC = — Jcos(#j — 9j — Aij), where 
6i is the phase of the superconducting order parameter at site i, J is the isotropic coupling 
energy between nearest-neighbor sites, and = (27r/0o) Si A • dl accounts for any applied 
magnetic field B through a vector potential A. The choice of triangular grid minimizes 



the periodic vortex pinning potential which is an artifact of any discrete lattice ||13|| . We 
evaluate the equilibrium and dynamical properties of Ji using respectively a standard Monte 
Carlo (MC) procedure fL4| , and solving the coupled Josephson equations with Langevin 



dynamics fl5f . For each plaquette p, we determine a vorticity vector n a (p) (a = x,y,z) 



from Y^ p mod[9i — 6 S — Ay, 2tt] = 2n[n a {p) — f p ]. Here the summation runs along the bonds 
belonging to the plaquette labeled p; and f p = Y^ p A^ / (2tt) . 
Previous calculations j|] have established that there is a first-order phase transition at 
high fields in the c direction (/ = 0/0 o = 1/6, where is the flux per triangular plaquette, 
and 0o — hc/2e). At low fields (f = 1/25) on a simple cubic lattice [H|, there is an apparent 



double transition: the first stage is a melting transition, followed at higher temperatures by 
a loss of phase coherence parallel to the c axis (as measured by j Z z)- By contrast, Jagla 
et al found that for / = 1/6, the transitions are separate only in the presence of quenched 
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We first show that a similar double transition occurs in the stacked triangular lattice. 
Fig. ^ shows j zz (T), the c-axis resistance R C (T), and the first order Bragg peak intensity 
/(Gi,T) (inset) as a function of temperature T for / = 1/24. Gi is taken as one of the 
six equivalent shortest reciprocal lattice vectors in the ab plane. We use periodic boundary 
conditions in all three directions and a lattice which is a parallelopiped of size 24 x 24 x 24, 
i. e., 48 vortex lines. Clearly, I(G\) and 7 ZZ go to zero at quite different temperatures, 
suggesting a double transition. The lower temperature, denoted T m ~ 1.55J, corresponds 
to the loss of solid-like structural order. The upper transition occurs at T#, ~ 2.0J. The 



behavior T m < T < Ti is characterized by logarithmically slow relaxation of ^y zz to its non- 
zero equilibrium values over (5 — 10) x 10 5 MC sweeps. The small high-T tail of ^ zz depends 
on sample thickness, and the data shown in the figure represent asymptotic behavior up 
to N z = 48. Starting from typical equilibrium configurations thus obtained, we carried 
out dynamical calculations within the resistively shunted Josephson junction (RSJJ) model 
|l5|Jl6|| . Using periodic boundary conditions, we injected a uniform current (0.083/ c per 
grain) through the lowermost xy-plane and extracted through the xy-plane in the middle 
of the sample, effectively creating a periodic pattern of alternating current flow along the 
z-axis. The resulting dissipation R c is sensitive mainly to the presence of net transverse 
vorticity, and its fluctuations. We also calculated R a b, analogous to the inplane resistivity, 
again using a similar periodic arrangement. This exerts a shearing force on two halves of the 
lattice and probe directly the shear rigidity of the lattice on length scales of our simulation 
box. The qualitative behavior of R a b is very similar to R c , i. e. it displays a sharp increase 
atT, [TJ. 

Fig. |2| shows the top-to-bottom vorticity density- density correlation function gtb{x, y) ~< 
n z (x,y,z = N z /2)n z (0, 0, 0) > as a function of T for / = 1/24. Note the periodic lattice 
of peaks at T/J = 1.5, the presence of only a strong central peak at T/J = 1.8, and the 
absence of peaks at T/J = 2.2. The strong central peak suggests that at T/J = 1.8 the 
vortex lattice has melted into a line liquid with a finite line correlation length along the field 
direction. This separation between the melting and ^ zz transition is a low-field phenomenon: 
at / = 1/6, the I(Gi) and j zz vanish at nearly the same T within our resolution. 

To probe the vortex density change associated with the transitions, we must abandon 
periodic boundary conditions, because these enforce a fixed density in the system, whose 
value is dictated by the applied field. We therefore adopt free transverse boundary conditions 
fi~8|| , while still maintaining periodic boundary conditions in the z direction. This procedure 
allows the outermost surface of the vortex ensemble to have a variable position. 

We first discuss our results in 2D, considering / = 1/24 on a triangular grid of different 
sizes with N = 26, 52, 78 and 100. One might think of taking the "magnetization" M z as 
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the average net vortex density n = / n z (r)dr/A, where n z (r) is the local vortex density and 
A is the total area. However, M z defined in this way suffers from spurious boundary effects, 
arising from the depletion of vortices near the boundaries in the lattice phase. This surface 
effect |T!J vanishes for large samples as 1/ \J~A~. Instead, we look at the local Voronoi cell area 
Ai, i. e., the area of the generalized Wigner-Seitz cell for vortex i ||. A local vortex density 
at a point R in the simulation box may then be defined as n(R) = J2i ^(R- ^ -A\)/-A^i where 
<5(R G Ai) gives 1 if the point R lies in the Voronoi cell associated with vortex i. Next, 
the local magnetization M z which we equate to the bulk average density < n > is calculated 
from < n >= X^ec x as the average of inverse Voronoi area for vortices lying within a 
measurement area C suitably distant from the sample boundary. 

For T > 0.5 J, there are, besides the field-induced vortices, numerous thermally induced 
vortex-antivortex dipole pairs, which must be eliminated before this procedure is applied. To 
do this, we pair each antivortex with the nearest vortex, and remove them from the count. 
Since most such dipole pairs are much smaller than 1/ \J< n >, this criterion is justified. 

For / = 1/24 the 2D melting transition occurs near = 0.09 J, as determined by the 
vanishing of the Bragg peak /(Gi). The zero-field XY (i. e., in 2D, the Kosterlitz-Thouless) 
transition occurs at about T^ T = 1.6J. Superconductivity in 2D is destroyed at T%f. As 
shown in Fig. || < n > /n (where n = y/3f/A is the nominal density) displays a sizeable 
increase across T^ T . This increase is independent of system size(checked up to 100 x 100 
grains). As further evidence that the boundaries are not producing this increase, we note 
that we obtain the same results for < n > for different choices of measurement areas, C\ 
and Ci- 

The upper panel of Fig. 3 shows that the local density fluctuates progressively more as T 
increases: near T m , the rms width of the local bond length fluctuation is about 0.15, a value 
close to the Lindemann number. The increase in density at high temperatures arises mainly 
from an increase in the frequency of short bonds, which occur throughout the sample as is 
clear from the panel. We have also observed at T m a sharp jump in the density of topological 
lattice defects, such as disclinations. If these defects are predominantly of a particular type 
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such as lattice vacancies, we would expect a corresponding jump in vortex density. However, 
we do not observe a sizeable increase in M z at T m within our resolution of about 2 %. 

As the vortex liquid is further heated, the vortex density increases to a value well above 
that of the solid. This increase is due to the excitation of vortex-antivortex pairs charac- 
teristic of the / = Kosterlitz-Thouless transition. The excited pairs screen the repulsion 
between the field-induced vortices, lowering the effective chemical potential for field induced 
vortices and making the system more compressible. 

Fig. |] shows the corresponding behavior in 3D, at / = 1/24 and 1/6. Once again 
ti/tiq ~ 1 below the zero field 3D XY transition at Txy, increasing markedly near T/J ~ 2.0 
for / = 1/24 and near 1.15 J for / = 1/6. For / = 1/6, this increases occurs over a 
temperature range AT < 0.05J, just at the first order phase transition where both lattice 
order and phase coherence parallel to the field are destroyed (T m ~ T £ for / = 1/6). The 
relative change in vortex density, Sn/no, is less than 7 % for / = 1/6. The entropy as 
calculated from f T C/TdT displays a corresponding jump AS ~ 0.11&B per grain (or 0.3/cg 
per vortex) as shown in the upper panel. 

For / = 1/24, the change in n/riQ is spread over a wider temperature range (AT ~ 0.7 J), 
and the midpoint of the increase in u/uq falls at Tg_. The overall field change SB/B ~ 
15%. The entropy does not show a recognizable jump near T m . For this open-boundary 
system, T m < 1.35J, about 10% below the value T m ~ 1.55J found with periodic boundary 
conditions. At / = 1/6, the two transitions are indistinguishably close, and any possible 
change associated with melting is overshadowed by the larger increase in n due to screening 
by vortex loops at for this dense system. Thus, for / = 1/24, there is a phase between 
T m and Ti with no lattice order, but yet with a diamagnetic moment almost as large as in 
the vortex lattice phase. At high temperatures, the vortex density rises substantially above 
the solid phase. 

This overall increase in M z is consistent with experiments, but the rise is associated, not 
with melting itself, but rather with the disappearance of phase coherence in the z direction 
{izz 0), which is also accompanied by dramatic changes in the transport coefficients, R c 
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and R a b- In 3D, this increase is associated with the proliferation of vortex loops, caused by 
the 3D XY fluctuations. The loops screen the repulsion between vortex lines, allowing them 
to come closer together and the net vortex density to increase. In 3D, the screening also 
encourages collision and fusion of field induced vortex lines [2(J, causing them to lose the 
longitudinal phase coherence at T(. 

The rise in M z for / = 1/24 in our model would actually occur over a very narrow 
temperature range in the real superconductor. This is because J in YB^CusOr-^ is strongly 
temperature-dependent near T c with J oc 16 J^(r) 2 • Using the simple mean field form for the 
temperature dependence of A(T) and material parameters for YB^CusOy-^, we estimate 
AT = [T e - T m ] < 150mK for / = 1/24, and AT < lOmK for f = 1/6. Thus the rise in 
M 2 , in fact, could occur over a very narrow temperature range in high-T c materials, in good 
agreement with experiments . 

Finally, we consider whether the upper transition in 3D might be a true, second phase 
transition, and not just a crossover. Numerically, there are several results suggesting a 
separate phase transition at Tg_. First, ^ zz — > at T# > T m |14|,0,pTJ, and dissipation due to 



extensive flux line cutting sets in both parallel and perpendicular to the field at Tj > T m fll? 
An analysis of vortex configurations also reveals that infinite clusters of colliding vortex lines 
and loops first appear at T [0 reminiscent of a percolation transition. Recent numerical 
simulations using London vortex loop model suggest a similar phase transition, separate 
from the melting transition p2[23[ . 



Tesanovic has considered the low field regime of the mixed state Jj]]. By a singular gauge 
transformation introducing fictitious antivortices, the original problem is transformed into 
one in which the new order parameter experiences no applied field but a fluctuating gauge 
field. The XY-type transition of this new order parameter(the so-called "$ transition") is a 
possible candidate for the transition at Tg for / < 1/24. We emphasize that increase in M z 
at Te results from the onset of screening by the background vortex loop plasma. It appears 
that a similar transition between a superfluid and metallic state exists in a 2d quantum 



bosons interacting with strong transverse gauge fluctuations p4| , |23 



In conclusion, we have calculated the change in magnetization and other properties of 
the flux lattice in a model high-T c superconductor near and above the melting transition at 
relatively low fields ||25|| . The transition proceeds in two stages, a lower melting transition 
being followed by a broader transition or crossover associated with proliferation of vortex 



loops in 3D, or vortex-antivortex pairs in 2D [ ZE| . The magnetization increases sharply at 
this second stage, in good agreement with recent experimental data. 

This work was supported by DOE Grant DE-FG02-90 ER45427 through the MISCON, 
and by NSF Grant DMR94-02131. 
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FIGURES 




T/J 



FIG. 1. Helicity modulus ^y zz (Filled circles); c-axis resistance R c (open circles) in units of 1/to 
where to = h/(2eJR). The error bars for R c show the standard deviation when the averaging 
interval was varied from 200io to 600to in steps of 10to- The inset shows the integrated first-order 
Bragg intensity 1{G\) of the field induced vortices, for / = 1/24. Error bars represent rms deviation 
over 7 different configurations 7000 MC steps apart. Lines are guides for the eyes. 
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FIG. 2. Top-to-bottom vortex density-density correlation function gtb{x,y) for T/J = 1.5,1.8 
and 2.2 in a 24 x 24 x 24 system with / = 1/24. T m ~ 1.55 J and T e ~ 2.0 J. The lower panel shows 
typical vortex configurations. 
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FIG. 3. Upper paneklocal density distribution n(R), plotted as a function of distance R from 
the center of the plane for T/J = 0.08, 1.05, and 2.05. Lower panel: normalized bulk vortex density 
vs. T for three different sample sizes at / = 1/24. Right panel: results of Delaunay triangulation 
for field induced vortices at T/J = 1.05. For results in left panel, averages were taken over vortices 
within two different bounded areas, C\ and C<i- 
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FIG. 4. Entropy S (upper panel) and normalized vortex density (lower panel) for / = 1/24 
and 1/6 in 24 x 24 x 12 system with open boundary conditions. The error bars denote the rms 
deviations from layer to layer. 
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